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We study the effect of dark matter (DM) particles in the Sun, focusing in particular on the 
possible reduction of the solar neutrinos flux due to the energy carried away by DM particles from 
the innermost regions of the Sun, and to the consequent reduction of the temperature of the solar 
core. We find that in the very low-mass range between 4 and 10 GeV, recently advocated to explain 
the findings of the DAMA and CoGent experiments, the effects on neutrino fluxes are detectable 
only for DM models with very small, or vanishing, self-annihilation cross section, such as the so- 
called asymmetric DM models, and we study the combination of DM masses and Spin Dependent 
cross sections which can be excluded with current solar neutrino data. Finally, we revisit the recent 
claim that DM models with large self-interacting cross sections can lead to a modification of the 
position of the convective zone, alleviating or solving the solar composition problem. We show that 
when the 'geometric' upper limit on the capture rate is correctly taken into account, the effects of 
DM are reduced by orders of magnitude, and the position of the convective zone remains unchanged. 
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I. INTRODUCTION 



The identification of Dark Matter (DM) is one of the 
most important open problems of modern cosmology and 
particle physics. Among the many DM candidates pro- 
posed in the literature, Weakly Interacting Massive Par- 
ticles (WIMPs) are the most popular, as they naturally 
achieve the appropriate relic density and naturally arise 
in extensions of the Standard Model of particle physics 
that will soon be probed in accelerators. WIMPs are cur- 
rently searched for in a variety of experiments (see Ref. |l| 
for recent reviews on particle DM, including a discussion 
of ongoing direct, indirect and accelerator searches). 

Here, we focus on the effect of WIMPs on stars, and 
more specifically on the Sun. In general, WIMPs cross- 
ing a star can in principle scatter off baryons to veloci- 
ties lower than the star escape velocity and be therefore 
gravitationally trapped by the celestial body, the effi- 
ciency of the process depending on the WIMPs scatter- 
ing cross section and WIMPs ambient density. Once cap- 
tured, WIMPs can still scatter with the star nuclei and 
therefore transfer energy inside the object. In addition, 
WIMPs can eventually annihilate, providing therefore an 
exotic source of energy. These effects have been investi- 
gated in the past, in particular in the context of WIMPs 
'cosmions' in order to solve the solar neutrino problem 
@, y] • This solution has then been discarded in favour 
of neutrino oscillations and the combination of the cos- 
mion scattering cross-sections and masses suggested by 
the solar neutrino problem is nowadays excluded by di- 
rect detection constraints. 

Still, the intriguing possibility that WIMPs modify the 
properties of the stars has recently prompted new interest 
on the subject. Recent works have focused in particular 



on compact objects |J, Q, main sequence stars at the 
Galactic center [fj, @] and first stars pfl-KL7| . In all these 
studies the celestial objects are placed in environments 
with high DM density, thus enhancing the numbers of 
WIMPs captured by the star and therefore the WIMPs 
transport and annihilation energy budgets. However, in 
these scenarios it is not possible to set robust constraints 
on the DM cross sections and mass because of the large 
uncertainties on the DM densities and/or the lack of pre- 
cise enough observations of the star targets. 

The situation is quite different for the Sun since its 
properties have been measured with good precision. Be- 
cause of that, the Sun can be used as a diagnostic tool 
to test small modifications of its structure induced by 
DM. The most recent works in this direction have fo- 
cused on modifications of the solar neutrino fluxes and 
helioseismology data [18h2C)( . Here we perform a com- 
plete and self-consistent calculation of the Sun evolution 
inside the galactic DM halo. In particular we extend 
previous works considering light DM candidates, with 
masses in the ran ge s uggested by DAMA [21| CoGent 
|2a | and CDMSII [23J experiments. We consider both 
standard annihilating WIMPs models and scenarios with 
negligible DM annihilations. We then move on Self In- 
teracting Dark Matter (SIDM) models [Hj], which has 
recently been invoked as a solution of the solar compo- 
sition problem [2[|. Finally, we focus on Inelastic DM 
models |26j . 

The paper is organised as follows: in Section |H] we in- 
troduce the DM models used in our calculations. In Sec- 
tion |TlT] we present the GENEVA stellar evolution code. 
In Sec II VI we discuss how the solar neutrino fluxes can be 
used to constrain the DM parameter space. We present 
our results for standard annihilating WIMPs in Section 
IVl for asymmetric DM models in Section IVII for SIDM 



models in Section IVIII and for Inelastic DM in Section 
IVIIII Finally we summarize our conclusions in Section 
IIXI Technical details and the formulas implemented in 
the code are in Appendix |A"1 



II. THEORETICAL SETUP AND NUMERICAL 
CODE 

Galactic WIMPs inevitably stream through any celes- 
tial object; those which scatter with a nucleus of the 
celestial object (in the following we will limit ourselves 
to consider stars) loose energy in such scatter, and if the 
energy loss brings the velocity of the WIMP below the 
escape velocity, the WIMP is "captured", i.e. it becomes 
bound to the star. The formalism to compute the rate 
at which WIMPs are captured by a star has been object 
of extensive studies in the 80's, and here we adopt the 
established results from 271. The evolution of the total 



number of WIMP particles N x inside the Sun (or any 
celestial object) can be written as 



N X = C 



2ANt - EN^ 



(1) 



where C is the particle capture rate over the Sun, A is the 
annihilation rate and E the evaporation one. We discuss 
the details of the time evolution of the WIMP number 
density in the Appendix \K\ and we limit ourselves here to 
discuss the parameters entering in the calculation and the 
implementation of the theoretical setup in the GENEVA 
numerical code for stellar evolution. 



A. WIMP parameters 

a. DM density. The number of WIMPs captured by 
a star depends on the WIMP scattering cross section off 
nuclei, ambient density and velocity distribution (see the 
Appendix for a detailed discussion). The DM density 
distribution in the inner regions of the Galaxy is poorly 
constrained. Semi-analytical studies and simple extrapo- 
lations of the results of N-body simulations predict high 
DM densities at the center of galactic halos, therefore 
suggesting stars at the galactic center [y, 0] and at the 
center or primordial halos [9rfl7| as promising targets to 
detect the effects of WIMPs on stars. 

In this paper we focus on the Sun, which being ~ 8 
kpc distant from the galactic center, is expected to be 
placed in a much lower DM density environment than 
the galactic center. Recent studies J28| - [3C| have shown 
that current observations constrain the local DM den- 
sity, p x , within a factor two, assuming spherical DM pro- 
files, while significant deviations may occur for oblate 
DM distributions. In the following, we adopt p x = 0.38 
GeV cm -3 , the value obtained in [28( for the benchmark 
Navarro-Frenk- White density profile. This value was ob- 
tained under the assumption that the DM profile is spher- 
ically symmetric, while recent numerical simulations that 



take into account the effect of baryons on the DM density 
profile suggest that the true DM density is at least ~ 20% 
higher, and that it is affected by large systematic uncer- 
tainties due to the unknown position of the Sun in the 
triaxial DM potential of the Milky Way J3l|. Further- 
more, the presence of a DM disc, corotating with stars 
and henceforth "colder" with respect to the solar motion, 
could enhance the capture rate in the Sun up to a fac- 
tor ten [32|] ■ Our conclusions can be trivially rescalcd for 
these cases, which however we do not take into account 
for sake of conservativeness. 

b. WIMP-nucleon scattering cross section. The 
WIMP scattering cross section off nuclei is constrained 
by direct detection experiments. The strongest bounds 
on Spin Indipendent (SI) interactions are currently set 
by CDMS-II HI, XENON100 HI and CoGeNT [H 



which impose crsi < 3 10 



-40 



4 10 
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cm in the 



mass range m x ^5—10 GeV. The strongest limit on 



Spin Dependent (SD) interactions is from PICASSO [34|: 
o~sd < 4 10~ 36 — 4 10~ 37 cm 2 in the same mass range. 
An evidence for an annual modulation signal has been 
reported by the DAMA collaboration and subsequently 
confirmed by the upgraded apparatus [2ll . |35| . |36( . A sim- 
ple interpretation of this signal in terms of WIMP clastic 
scattering off nuclei of the target material can hardly be 
reconciled with the constraints from the other DM detec- 
tion experiments. In this framework, in fact, the bulk of 



the region of the a - 



parameter space favored by the 



DAMA signal is excluded by the null results of the other 
experiments, both considering SI and SD interactions. 
Still, there remain small surviving regions of the pa- 



rameter space for low WIMP masses, 



10 GeV. 



We address the reader to specific studies for more details 
(e.g. [23, Hal an d references therein). Recently, two ex- 
periments, CDMSII H| and CoGeNT [H[, have reported 
an excess in their recoil energy spectra, even if the sta- 
tistical significance of these signals is very weak. Inter- 
estingly, if interpreted in terms of elastic WIMPs scat- 
tering, these results would point to light WIMPs, with 
masses and scattering cross-sections close to the DAMA 
region [38|, [39| . This scenario has been challenged by the 
very recent upper limits on SI interactions obtained by 
XENON100 [33(. However, there are subtle experimen- 
tal effects that may cast some doubts on the low-mass 
end of this result, and as we write, the debate between 
the CoGent and Xenon collaborations is still open. Since 
the WIMP energy transport inside the star is indeed en- 
hanced for light WIMPs, in Sec lVIl we investigate the 
effects of these particles captured by the Sun on its struc- 
ture and obscrvables. 



B. Dark matter models 

c. Standard WIMPs. Standard WIMPs achieve the 
DM abundance inferred from cosmological observations 
for (crv) ~3x 10 -26 cm 3 s _1 . Although significant devia- 
tions from this value are expected for example in presence 



of Sommcrfcld enhancements, efficient coannihilations or 
for non standard cosmologies (see Ref. [l[ for more de- 
tails and references), deviations of few orders of magni- 
tude from the value quoted above are not relevant for the 
effects of the WIMPs in the Sun, as we shall see in Sec. 
|A"1 Indeed, the WIMPs annihilation and energy trans- 
port are independent from (av) once an equilibrium be- 
tween scatterings and annihilations is reached. This oc- 
curs within a timescale larger than the age of the Sun only 
for (av) < 10 " 



-33 



cm 3 s _1 , after Equation IA7I Therefore, 
for the rest of the paper we fix (av) = 3 x 10~ 26 cm 3 s _1 
when we consider annihilating WIMPs models. 

d. Asymmetric dark matter DM annihilations can 
be neglected cither in models with very small or van- 
ishing annihilation cross sections or in presence of an 
asymmetry in the DM sector between particles and anti- 
particles. This possibility, realized in the so-called asym- 
metric DM models, is particularly interesting since these 
DM candidates can in principle have weak-scale interac- 
tions and therefore sizable scattering cross sections off 
baryons, despite they do not annihilate. Concrete real- 
izations of this idea are for example models where the 
dark sector contains a conserved U(l)x symmetry, ana- 
logue of the baryon number, responsible for the stability 
of the lightest particle in the DM sector. If this quantum 
charge is shared between baryons and DM it can link the 
asymmetries in the two sectors and this may naturally 
explain why the baryons and DM abundances are of the 
same order of magnitude (see e.g. |40H44| ). As noticed 
above, all DM models for which, the value of the annihi- 
lation cross section implies an equilibrium time r© > £© , 
namely (av) ;$ 10~ 33 cm 3 s _1 , can be considered asym- 
metric from the point of view of WIMPs accretion. 

e. Self- Interacting dark matter. In the scenarios de- 
scribed before, the capture of WIMPs in the star is in- 
duced by the interactions of DM with the nuclei. In ad- 
dition to that, also the elastic scattering of star crossing 
WIMPs off previously captured WIMPs can contribute 
to the DM capture rate. For WIMPs self-interactions 
of the order of the WIMPs scattering cross sections of 
baryons, like in the models previously considered, this 
contribution is completely negligible since the number of 
target nuclei in the Sun is much higher than the number 
of trapped WIMPs. However, this is not longer true for 
enhanced self interactions, as predicted in Self Interacting 
Dark Matter (SIDM) models. In Ref. [HI it has indeed 
been shown that for extreme values of self interaction 
cross sections and for (av) < 10 -27 cm 3 s _1 , the annihi- 
lation rate can indeed be boosted by a factor of order 
0(100), consequently enhancing the prospects for detec- 
tions of high energy neutrinos produced by WIMPs anni- 
hilations in the Sun. SIDM models have been proposed 
in Ref. [24j to solve some discrepancies of the collision- 
less Cold Dark Matter (CDM) model, in particular the 
excess of substructures found in N-body simulations of 
CDM halos with respect to observations and the conflict 
between the predicted cuspy CDM density profiles and 
the observed cored halos of LSB galaxies. Concrete par- 



ticle physics realizations of SIDM models includes theo- 
ries with strongly interacting dark sector W2M4J] , Q-balls 
j24l l46j , quark-gluino bound states [43, |48( and mirror 
DM models |49l - (52| . Observations of galactic and galaxy 
clusters halos constraint the size of DM self interactions. 
Here we report the bound obtained from the bullet clus- 
ter as a robust upper limit on the SIDM self interaction 
cross sections [53( : oy v /m v < 2 x 10~ 24 cm 2 GcV _1 . We 
address the reader to Ref. [54[ for an updated discussion 
on the constraints on SIDM. 

In the following, since the main focus of the paper is 
the effect of WIMPs on the Sun, all stellar quantities we 
show are snapshots taken at the age of the Sun. 



III. THE NUMERICAL CODE 

In order to perform an accurate analysis of the prob- 
lem we have implemented the capture, annihilation and 
energy transport of WIMPs in the GENEVA code (for a 
description of this code, see |55J|). 

The implementation was done in the following way: 
from the structure of the stellar model, we can compute 
at each time step the number of WIMPs in the stellar 
model, the energy released by annihilation and the en- 
ergy transferred by the WIMPs. At their turn, these 
energies modify the structure of the model. Thus the 
problem has to be resolved self-consistently. To do that, 
the rate of WIMPs annihilation and transferred energies 
are added in the equation describing the conservation of 
the radiative energy. The code then, through an iterative 
procedure, looks for the values of the pressure, temper- 
ature, radius and luminosity at different depths in the 
stellar model which match for the WIMP modified stel- 
lar structure equations (usually a stellar model consists 
of 900 layers from the surface to the centre) . 

With this numerical tool, we computed the evolution 
of one solar mass stellar models from the Zero Age Main- 

9 

Sequence (ZAMS) until the solar age (4.57 xlO years) 
for various prescriptions for the WIMPs energy transport 
and annihilation. As initial composition we chose the 
following values: mass fraction of hydrogen, X=0.72, of 
helium, Y=0.266, and of the heavy elements, Z=0.014. 
The distribution of the heavy elements was taken as given 
by [56[. We did not include the effects of microscopic 
diffusion in these models, although this effect would be 
required for the computation of tailored models for the 
Sun [e.g. [53]. The inclusion of this effect would however 
not change the conclusions of the present work, which 
are deduced from detailed comparisons of models with 
exactly the same physics except for the inclusion or not 
of the WIMPs. Since, microscopic diffusion would have 
to be accounted in both models, it would not change in 
a significant way the relative effects obtained here. 



IV. DIAGNOSTIC TOOLS 

The modification of the stellar structure produced by 
the WIMPs induces changes in the frequencies of stellar 
oscillations modes and in the neutrino fluxes. The first 
signature can in principle be observed with helioseismic 
measurements [18|-|20| , however, the neutrino flux is much 
more sensitive to the variation of temperature and den- 
sity profile of the innermost regions of the Sun, and hence 
it is a much more powerful diagnostic tool. 

For the standard solar model, Ref. |58| indicates that 
the 8 B neutrino flux varies as T 25 , therefore a 1% tem- 
perature change will produce a 25% change in the 8 B 
neutrino flux. However, as it has been in shown in Ref. 
J2G | , this simple scaling law is not valid to describe the 
peculiar modifications of the temperature profile induced 
by the WIMPs. Therefore, the use of a stellar code is 
mandatory to study the effects on DM on the solar neu- 
trino fluxes. 

The distribution of WIMPs inside the Sun is crucial to 
determine the modifications of the neutrino fluxes. Once 
captured, WIMPs get redistributed inside a small spatial 
scale, r x , of the order of I0 9 cm - 0.01 R Q , for a WIMP 

mass of 100 GeV (r x scales as m x , see Appendix lA2p . 
The WIMPs energy transport will be therefore more ef- 
ficient in the innermost regions of the Sun core, as it can 
be appreciated in the left panel of Fig. [2] (see next sec- 
tions for more details.) The 8 B and 7 Bc neutrinos, which 
are mostly produced at ~ 0.04i? Q and ~ 0.06i? Q , will be 
more affected by the presence of WIMPs than the pp neu- 
trinos. In fact, although the pp neutrinos are the most 
abundant solar neutrinos, they are mainly produced at 
~ O.li?©, thus is a region well outside the one affected 
by the WIMP energy transfer. For the same reason and 
considering the experimental uncertainties in the deter- 
mination of 8 B and Be neutrino fluxes, we conclude that 
the 8 B neutrino flux is the best diagnostic tool in order 
to test the effects of WIMPs on the Sun. 

The 8 B flux has been determined with good accuracy 
bySNO ||: 

<f> B = 5.046i£^(stat)+oi23(syst) 10 6 cnrV 1 . 

For completeness, we report also the 7 Be neutrino flux 
inferred by Borexino[60j: 

<f> Be = (5.18 ±0.51) x 10 9 cnrV 1 . 

The solar model we use, as described in IIII1 predicts at 
the solar age t Q : 



4.56 x 10 b em's" 



= 4.47 x 10 9 cm _2 s 



values that are well in agreement with the experimen- 
tal results to within the theoretical uncertainties of solar 
model calculations |6lTl63| . 

Despite its great success in explaining a large variety 
of observations, the standard solar model suffers nowa- 
days from the so-called solar composition problems. Re- 
cent analysis points to a lower surface heavy element con- 
tent than previously thought (see [64| for recent results) 



and solar models incorporating these revised metallic- 
ities conflict dramatically with helioseismological mea- 
surements, in particular right below the solar convective 
envelope. With the solar abundances of [64j, the radius 
of the convective zone, Rcz, is more than 10 a higher 
than the measured value: Rcz = 0.713 ± O.OOl-R© (see 
Ref. [6 If for an analysis of comparison of different solar 
models with hclioseismology measurements.) High and 
low heavy elements models also produce differences on 
4> B and 4> Be respectively of ~ 20% and 10%. However, 
the theoretical uncertainties obtained within a partic- 
ular solar model (with high or low metallicity estima- 
tions) are sig nificantly smaller, ~ 10% and 6% for 
and <b 



02. (k 



Considering the theoretical and ex- 
perimental uncertainties on <j> B , we study the region of 
the DM parameter space where significant deviations on 
(p B are found (the maximum allowed deviation of the 8 B 
flux depends on the value of the theoretical uncertain- 
ties considered. We define this thereshold in Sec lVI[) . As 
reference value, we consider the neutrino flux obtained 
from our solar model without DM particles and in the 
Sec. I VII we define the maximum deviations compatible 
with present data. 

We now discuss the results for standard WIMPs, and 
provide a description of the role of WIMP energy trans- 
port and annihilation. We present separately in the next 
two sections the results for two special classes of models: 
asymmetric and self-interacting DM. 



V. RESULTS FOR STANDARD WIMPS 

Since the trapped WIMPs are confined in small spa- 
tial scale inside the Sun, the annihilation of WIMPs in 
the star acts as a point-like source of energy, whose ef- 
fects have been studied in details in the context of small, 
and high mass stars, and for various initial metallicitics 
(d,0,la-[Il,|6i|). When the energy injected through DM 
annihilation is comparable with the nuclear one, a star 
finds its equilibrium stage at lower temperatures, thus 
stopping in a colder region of the Hcrtzsprung-Russel di- 
agram and prolonging its lifetime. Although the effect 
is in principle interesting, the required DM densities and 
cross-sections are such that the effects of an additional 
energy source due to the DM are totally negligible on the 
Sun, and anywhere else in the Galaxy but in its inner- 
most regions 0, [7J . The WIMPs in the star are confined 
into closed orbits and they inevitably scatter off the stel- 
lar material therefore redistributing the energy inside the 
star. The WIMPs mean free path is l(r) = 1/ ^~\ o"jrij(r), 
where the sum is performed over all the target nuclei, Oi 
is the WIMP scattering cross section off a certain nu- 
cleus and rii(r) is the nucleus number density profile in 
the star. For reference, the WIMP mean free path is 



10 



-.'56 



cm 



typically of the order ~10 cm for asD 
Once compared with the typical radius of the WIMP 
cloud r x =10 9 cm, this indicates that WIMP transport 
effects are non local. In our computations we span many 
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Figure 1: Comparison of the energy transported (etrans, blue lines) and injected (red) per unit mass and unit time by WIMPs in 
the Sun, with the nuclear e nuc i and gravitational e gra v energy density of the Star, for m x = 10 GeV (left panel) and m x = 100 
GeV (right). Here we plot the absolute value of e tT ans, the sign being negative at radii smaller than the dips of the blue curves 
at ~ 0.04 Rq (~ 0.014 Rq) in the left (right) panel, and positive otherwise. The sign of e grav is positive, this meaning a 
contraction of the star. 



orders of magnitude on the WIMP scattering cross sec- 
tions so we investigate both local and non-local regimes. 
We address the reader to Ref. Q, and references therein 
for further details. The equations implemented in the 
GENEVA code are summarized in the Appendix [A] for 
sake of completeness. 

In the non local regime regime, WIMPs scatters can ef- 
ficiently transport the heat from the inside of the stellar 
out to colder regions, thus operating toward a flatten- 
ing of the temperature profile, and a consequent read- 
justement of the entire stellar structure. Raising the 
WIMP-proton scattering cross section a , the capture 
rate (modulo the caveat on saturation of the optical limit, 
see the Appendix IA ip and the WIMPs scattering rate 
grows, therefore increasing the energy transported by the 
WIMPs, e t rans- The dip in the behavior of etrans in Fig- 
ure Q] corresponds to the change in sign of etrans'- at small 
radii WIMPs are "colder" than the baryons, and acquire 
energy with their scatters, which they deposit outward. 
Since the mean free path of WIMPs scales like I ~ l/oo, 
at large scattering cross sections the WIMPs remain pro- 
gressively "trapped" in the interior of the star. This im- 
plies that the heat transport from DM particles becomes 
local and therefore e tra ns dramatically reduces. 

The WIMP mass also plays a role in determining the 
effects on DM in stars, as it can be appreciated in Equa- 
tions in Appendix [A] Lowering the DM mass goes in the 
direction of maximizing the transport effects, but also 
the evaporation rate. Above m x = 5 GeV, evaporation 
can be safely neglected and the mass of the particle acts 
to modify the radius and normalization of DM inside the 
star. 

In Figure [JJ we compare the rate of energy released 



per unit of baryonic mass by DM annihilation e an n and 
the rate of energy absorbed or released by DM trans- 
port etrans-, with the rate of energy produced by nu- 
clear reactions. The solid lines represent two bench- 
mark WIMP scenarios for the Sun: a local DM den- 
sity of p x ~0.38 GeV/cm 3 , a Spin-Dependent scattering 
cross section os£)=10 _38 cm 2 and WIMP mass and self- 
annihilation cross-section compatible with a DM ther- 
mal production scenario. Since in these computations 
the WIMP annihilations and transport energies are neg- 
ligible with respect to the energy provided by nuclear 
reactions, the structure of the Sun is not affected by the 
presence of the WIMPs. As a consequence of that, the 
solar neutrino fluxes are unchanged with respect to their 
values in our standard solar model without DM parti- 
cles. As stated before, the WIMP energy transport is 
more efficient for light WIMPs: in Figure [TJ it dominates 
over the annihilation energy for a 10 GeV WIMP while 
the opposite trend occurs for a 100 GeV WIMP, at least 
in the inner region of the Sun. However, due to the dif- 
ferent scaling of e an n and etrans with p x , at high DM 
densities the annihilation energy is always the most rel- 
evant source of energy produced by the WIMPs in the 
star. An example of that is given by the dashed curve in 
Figure [TJ where we raise the DM density of two orders of 
magnitude, p x =38 GeV/cm 3 . 

Rescaling our curves in Fig. [TJ using Eqs. [TJ IA3I and 
IA61 it can be noticed that energy injection by DM an- 
nihilation is comparable with the nuclear energy for DM 
densities of order p x ~10 GeV/cm . This value is indeed 
the typical DM density for which the behaviour of main 
sequence stars starts to be severely modified by WIMPs 
annihilations [fj, 0] • 
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Figure 2: Impact of asymmetric DM on the Sun. Left Panel: Temperature profile of the Sun for different SD scattering 
cross-sections at £0 = 4.57 Gyrs. Central Panel: Density profile (baryons only) inside the Sun. Right Panel: Energy transport 
(meaning of the curves as in Fig[l}. The sign of e grav is positive. 



In conclusion, the Sun's global properties are not dra- 
matically changed by self-annihilating WIMPs, and no 
solar observable can be used as a good diagnostic for 
these models. As we explain in the following Sections, 
things arc quite different for other DM scenarios. 



VI. RESULTS FOR ASYMMETRIC DM 

If WIMPs are asymmetric, Equation [T] reads: 



N, 



C-EN X , 



(2) 



and particle will continue accumulating in the center of 
the star, since their abundance will not be limited by 
annihilation. 

As before, DM particles rapidly thermalize, therefore 
the shape of the WIMPs distribution stays as in the an- 
nihilating case whereas the normalization N x will result 
modified. The number of WIMPs inside the Sun, ne- 
glecting evaporation, will simply be N x = Ct Q , with 
£0 the age of the Sun. In the annihilating case, af- 
ter a transient r x , an equilibrium between annihila- 
tions and capture is reached and N x stays N x — Ct x . 
Comparing £q = 4.57 x 10 9 yr with the typical equi- 
librum timescale in a vanilla WIMPs scenario, r x ~ 
10 6 yr ((crt;)/lCr 26 cm^" 1 )" 1 / 2 for m x = 10 GeV, it is 
evident that the number of WIMPs trapped in the Sun 
is significantly bigger in asymmetric models. 

In order to study the effects of asymmetric DM, we 
have computed the variations of neutrino fluxes produced 
by the heat transport of WIMPs in the Sun. We evolved 
the Sun from the ZAMS up to its current age £q with 
the GENEVA stellar code. As previously discussed, the 
WIMPs effects are accounted for self-consistently at each 
time-step; if, for instance, energy is evacuated by the 
WIMPs in the very central layers, the structure will react 
by allowing the central layers to contract. The present 
numerical approach is therefore more accurate and self 



consistent than methods in which the effects of WIMPs 
are deduced from the structure of standard solar mod- 
els computed without WIMPs effects, this, for at least 
two reasons: first because in the present approach the 
structure of the star can readjust itself to any changes 
produced by the WIMPs and second because, these read- 
justments are accounted for during the whole previous 
nuclear evolution of the Sun. 

The first two panels of Figure [2] shows the temperature 
and density profiles of the Sun in presence of asymmet- 
ric DM. Note the change in temperature with increasing 
SD scattering cross section: the temperature decreases 
at the center and (although this is difficult to appreci- 
ate in the plot) slightly increases close to the external 
edge of the stellar core. The reason of the decrease in 
temperature can be understood in terms of energy trans- 
ported away from the solar core. In general, any amount 
of energy removed from the core leads to its contraction. 
Naively, one would expect a warming of the internal re- 
gions of the Sun since part of the energy extracted from 
the gravitational energy reservoir goes into internal en- 
ergy. However, the energy carried away by the WIMPs 
is well above the energy released by the core contraction, 
resulting therefore in a cooling of the central regions. 

The right panel in Figure [2] demonstrates in fact that 
the energy transported by asymmetric DM can become 
extremely high, and it can significantly modify the struc- 
ture of the star. We show the results for two WIMP 
masses: at fixed asD the normalization of etrans is higher 
for a 5GeV particle than for a 7GcV one, and the sign 
inversion point of the energy transport term shifts out- 
wards. The shift in the sign inversion point between dif- 
ferent curves for the same DM mass is due to the feedback 
induced by the modification of the density profile, which 
starts becoming sensible at <jsd «10 -36 . 

In principle one should study how the structural mod- 
ifications modify the details of neutrino oscillations in 
the Sun and consequently the inferred values of the neu- 
trino oscillations parameters, i.e. Atoq and sui(6!q). A 
sizeable effect is actually expected for density variations 
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Figure 3: Solar neutrino fluxes in presence of asymmetric DM in the Sun. Left Panel: Differential 8 B and 7 Be neutrino fluxes 
as a function of radius. The curves referring to our solar model are normalized to unity. Right panel: isocontours of 25% (red) 
and 5% (yellow) (f> v B deviations with respect our solar model prediction in the m x -asD plane. The red area show the region of 
the parameter space with (f> B modifications larger than 25% and therefore in tension with present (fi B data. The weakening of 
the DM transport effects at high cross sections is due to the transition to local transport regime. See text for details. 



of the order of 10% at ~ 0.04i? Q , which corresponds to 
the region of maximal 8 B production [661 ] . However, at 
that radii, we find a decrease of the density below 2% 
even for models producing modifications of the 4> B larger 
than 20-25%. This modification is comparable with the 
uncertainties on the density profile in the inner regions 
of the Sun, which is of the order of 1% for R < 0A5Rq 

These considerations arc valid also for the 7 Bc fluxes 
and therefore modifications of the neutrino oscillation 
probabilities are completely negligible for our purposes. 

Left panel of Figure [3] shows the differential 8 B and 
7 Be neutrino fluxes as a function of the stellar radius, 
in presence of different WIMP models. As expected, the 
reduction of the neutrino production, due to the cooling 
of the baryons inside the Sun, is more efficient at small 
radii, where the WIMPs are concentrated. As noticed in 
Sec. IIV| this leads to a larger modification on the total 



B neutrino flux, (j> B , which is the integral over the whole 
Sun of the corresponding differential quantity plotted in 
the left panel of Figure [H than those on the 7 Be neutrino 

flux, r Be - 

To study the impact of these structural variations on 
the solar neutrino flux, we have performed a systematic 
study of the DM parameter space, varying the WIMPs 
scattering cross section and mass. For SI interactions we 
find sensible variations of the neutrino fluxes only for very 
large values of a si, already severely excluded by direct 
detection experiments, therefore we focus for the rest of 
the section on SD interactions. 



In the right panel of Fig. [3j we show in the 



CSD 



plane the isocontours corresponding to (f> B variations of 
25% and 5% with respect to our solar model without 



WIMPs. For masses below 4 — 5 GeV the evaporation 
becomes important and the number of WIMPs inside the 
star is strongly suppressed. 

Above m x = 20 GeV the WIMP transport starts to 
become inefficient and even for high scattering cross sec- 
tion (still in the non-local transport regime) the WIMP 
energy transport is non negligible only at the very center 
of the star, providing a local dip of energy. However, for 
increasing WIMPs masses the existing constraints from 
direct detection experiments become more severe, and 
the region of the parameter space able to produce size- 
able modifications of cj) B is already excluded. Because 
of that, we do not explore that region any further. At 
very high scattering cross sections (i.e. >I0 _33 cm 2 ) the 
WIMPs heat transport becomes more and more localised 
and the modifications of (f> B tend to decrease, consistently 
with what is described in Sec|V] This is also the reason 
of the non specularity of the exclusion curves in Figure [3] 
the low and high cross-section regions are characterized 
by different physics (non-local vs local transport effects, 
respectively) . 

Combining experimental and conservative (20%) the- 
oretical uncertainties we derive that modifications of the 
4> B above ~ 30% are excluded at 95% CL. The maximum 
modifications that we obtained in our computations are 
slightly below this level: further increasing usd, prob- 
lems are encountered in solving the stellar structure at 
ages well below the solar one. We have not further inves- 
tigated if these difficulties are merely a numerical artifact 
or are instead related to the non existence of a solution 
of the stellar structure. However, we notice that once 
we obtain variations of <f> B of the order of 20%, further 
small changes of ctsd induce rapid modifications of <f> B . 



Because of that, the isocontours corresponding to 30% 
<f) v B variations should be closed to the ones corresponding 

to (5 



'B 



25% shown in Fig. El 



Considering a more optimistic value for the theoretical 
uncertainties on the <p B predictions, i.e. 10 %, the thresh- 
old for exclusion at 95% CL is lowered to 18% variations 
from our solar model prediction. Most of the SD cross- 
sections inside the 25% region in Fig|3] are excluded by 
the direct detection experiments constraints, apart from 
a small region at low masses which is somewhat in ten- 
sion with those bounds. A reduction of the theoretical 
and experimental uncertainties on (f> B in the next years 
may in principle improve the sensitivities on cfsd- We 
show however in Fig|3]that considering a 5 % modifica- 
tion of (j) B the region of the m x — <jsd parameter space 
which can be probed enlarges very little. Our fin ding s 
are qualitatively consistent with previous results |3j, l20j . 



VII. RESULTS FOR SELF-INTERACTING DM 

As we have seen, since DM annihilations limit the num- 
ber of trapped WIMPs in the Sun, asymmetric DM mod- 
els are promising scenarios to look for modifications of 
the Sun properties. This is particularly true for the case 
of SIDM models, as it has been noticed in Ref. [231 ■ I n 
presence of only DM self-interactions and DM scattering 
off baryons, the number of DM particle inside the Sun is 
given by: 



iV-, 



C + C XX N X 



(3) 



where C xx is the DM capture rate through DM self 
interactions. The solution of this equations reads: 



N x (t) 



C 



C v 



g^XX* 



1) 



(4) 



and reduces to N x (i) = Ct for negligible self interac- 
tions, i.e. the formula used in Scc lVIl It can be noticed 
therefore that the numbers of captured DM particles is 
exponentially cnchanced in asymmetric SIDM models. 
The self-capture rate C xx has been obtained in Ref 



C\ 



3 v xx v 2 {r) 
2 x m v v 



Erf (t?) 
V 



(5) 



where the average of <j>(r) = v 2 (r)/v 2 (R*) over the 
WIMPs distribution gives (<j> x ) ~ 5.1. i?» refers to the 
Sun radius and the definition of the other quantities can 
be found in Appendix. 

An upper bound to the total WIMPs self-interaction 
arises when the sum of the short range WIMPs self- 
interaction cross sections over all the WIMPs targets 
equals the surface of the WIMPs distribution: 



In the equation above, we can safely consider that all 
the WIMPs arc localized inside the thermal readius r x . 
Once the geometrical limit is saturated at a time t with a 
number of WIMPs N(t) inside the star, we simply replace 
the combination a xx N x in Eq|3]with cr e // iX . The WIMPs 
number density for t > t is therefore given by 



a eff,: 



T XX N X 



N x (t)=(C + C xx )(t-t)+N(t) 

with C xx obtained replacing a xx in Eq. [5] with <J e ff,x- 

We have studied the effect of a population of SIDM 
models in the Sun, implementing in the GENEVA code 
cquations[3H5] To exemplify our results we focus on a DM 
mass m x = 7 GeV and we consider agi = 2 x 10~ 41 cm 2 
and cfsd = 10~ 36 cm 2 , i.e. values close to the bounds 
set by direct detection constraints. In Fig. 2] we show 
the evolution in time of the number of captured WIMPs 
normalized to the total number of baryons, N x /Nq for 
the two models here considered and for a self-interaction 
cross sections which saturates the bound from the bul- 
let cluster mentioned above. In the case of a pure SI 
interaction it can be noticed the change between the ex- 
ponential rise regime and the linear one, occurring when 
the geometrical limit is saturated. For SD interactions 
the same happens but since the geometrical bound is ob- 
tained earlier in the evolution this modification is difficult 
to notice. In both calculations we have not found any sig- 
nificant deviation of the solar structure at the Sun age. 
For the case of SD interaction, the WIMPs transport en- 
ergy etrans is at most a factor 2 smaller than the nuclear 
energy. Small radiative opacity variations occur only at 
r < 0.05i? Q . These small effects does not lead to any 
significant variations of the sound speed profiles in the 
whole Sun , and so of Rcz ■ 

In Ref. [25| . adopting a polytropic model for the Sun 
structure and the linearized solar model of Ref. [68| 
the authors concluded that for suitable values of SIDM 
parameters the boundary of the solar convective zone 
Rcz is decreased. They claimed that these modifications 
are such that the solar composition problem mentioned 
above can be solved, i.e. the position of the convective 
zone and hclioscismology data can be reproduced. 

Our analysis does not confirm their conclusions. In or- 
der to maximize the WIMP energy transport, we have 
artificially turned off the geometrical bound in our code 
and evolved different solar models. We stress however 
that this scenario is unphysical. We have studied both 
models with pure SD and pure SI interactions, leaving 
the DM mass and cross section at the values fixed be- 
fore. We have then varied the DM self-interaction cross- 
section in order to maximize the number of SIDM in the 
Sun. In the extreme cases we have found a dramatical re- 
duction of the core temperature, which produces a signif- 
icant decrease of <f> B , in agreement with what is obtained 
in Sec. IVIl At the center of the star, the luminosity is 
reduced, due to the evacuation of energy produced by 
the WIMPs, and the radiative opacity is increased, as a 
result of the decrease of the temperature and the increase 
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Figure 4: Number of WIMPs inside the star normalized to 
the number of baryons as a function of the age in the asym- 
metric SIDM scenario. Solid and dashed lines refer to calcula- 
tions which respectively implement or neglect the geometrical 
bound on the self interaction. 



of the density. The external zones are not significantly 
affected by the WIMPs so we do not see any significant 
change in the position of the convective zone. We con- 
clude that a population of WIMPs, being strongly local- 
ized at the center of the Sun, can not affect its external 
shells without at the same time completely change the 
internal structure of the star. On the other hand, strong 
modifications of the central solar structure traduce in 
dramatic changes on the solar neutrino fluxes and helio- 
scismology g- modes, as also demonstrated in Ref. [69| . 
and in general of any observable sensitive to the physi- 
cal conditions in that inner regions. The results of Ref. 
[69] appear in good agreement with our own and confirm 
that the presence of WIMPs inside the Sun can not mod- 
ify the solar sound speed profile in such a way to restore 
the agreement with hclioseismological data. 



VIII. INELASTIC DM 

The possibility that DM particles could scatter inelas- 
tically up to an excited state, has been recenly proposed 
in the literature in the framework of theoretically justi- 
fied models [26|, [70l472l | , in order to explain the modula- 
tion observed in the DAMA signal. The existence of two 
states splitted by an energy 5 introduces a threshold for 
the WIMPs scattering since only WIMPs with enough 
kinetic energy, Ek can be excited to the heavier state: 



with tojv the mass of the nuclei. Since the energy thresh- 
old depends on the type of nucleus considered, in the in- 
elastic DM scenario the scattering rate can differ dramat- 
ically between different targets, opening the possibility to 
reconcile the DAMA sig nal with the results of the other 
experiments [38|, l73l - l76| . Recent analyses favor WIMP 
masses of the order m x ~ 10 — 50 GeV and small mass 
splitting 5 ~ 30—130 KcV. We address the reader to the 
original literature for more details. 

The WIMPs capture rate in the Sun for inelastic DM 
models tends to be reduced with respect to the clastic 
case by the presence of the inelastic barrier. However, 
since a significant fraction of the kinetic energy lost by 
the WIMPs goes into excitation, less energy is trans- 
ferred to the nuclei, thus reducing the nucleus form-factor 
suppression. Explicit calculations have shown than, for 
certain choices of the DM parameters, the capture rate 
can be enhanced up to an order of magnitude [77I . |78| . 
Another important effect of the existence of the energy 
threshold S for scatter is the modification of the internal 
distribution of DM particles inside the star: as a conse- 
quence of the reduced number of scatters, they can't shed 
away their angular momentum entirely, and their orbits 
remain larger than in a standard, elastic DM case (see 
Fig lain J77J). 

Here we aim to understand if inelastic DM models can 
affect the Sun's properties, and if the diagnostics we have 
discussed in previous Sections allows to probe this class of 
models. First, we notice that, once captured, the veloc- 
ity of DM particles will be smaller than the local escape 
velocity. Considering an escape velocity at the Sun core 
of v = 1300km s^ 1 and taking a WIMPs speed at infinity 
u = 220km s~ the maximum kinetic energy of a trapped 
WIMPs will be E 7 k nax ~ 96.5 KeV(m x /100 GeV). Look- 
ing at Eq. |6l it can be realised that once captured, 
WIMPs can scatter only with elements heavier than hy- 
drogen, even considering a small mass splitting 5 = 30 
kcV. Simple kinematic considerations show that after 
a few scatterings the WIMPs will not have enough ki- 
netic energy to overcome the inelastic threshold. For ex- 
ample, for the same S = 30 keV, a WIMP with mass 
m x = lOOGeV will scatter in average ~ 5 times off iron 
or off helium. Therefore, at a given time, only a small 
part of the WIMPs population inside the star can ef- 
fectively scatter and transport energy Focusing on the 
most optimistic asymmetric case, the number of WIMPs 
above the kinematic thereshold, N' at given time can 
be obtained from: 



N' 



C 



N' T~ 



(7) 



where Tj„ is the average time needed for a WIMP to 
fall below the inelastic thereshold. This quantity can be 
estimate as: 



E k >S[l 



M N 



(G) 



N s 



w 



(8) 



10 



where I is the WIMP mean free path defined in Sec. [Vj w 
is the WIMP velocity inside the star and N sc is the num- 
ber of WIMP scatterings. Here we take a mean baryon 
density p ~ 10 2 g cm" 3 and we consider only scattering 
off elements heavier than helium, which in total account 
for ~ 1% of the Sun mass. We obtain therefore: 



2.3 x 10 5 s 



N sc 1318km s" 1 lOOg cm" 3 10~ 36 cm 2 



30 



P 



&SD 



(9) 

At the age of the Sun, the number of WIMPs trans- 
porting energy inside the star will be Cn n , which is many 
orders of magnitude lower than in the elastic case, Ct & . 
Considering the capture rate obtained with the most 
favored values of inelastic DM parameters, from the re- 
sults obtained in Sec. IVII we conclude that inelastic DM 
models do not produce any detectable modifications of 
the Sun properties. 



IX. CONCLUSIONS 

By the use of a stellar evolution code, we have stud- 
ied the modifications on the Sun structure induced by 
DM particles captured by the Sun during its lifetime. 
For standard WIMPs an equilibrium between the capture 
and the annihilation rates is reached in short timescales. 
After this transient period, the number of WIMPs inside 
the Sun is constant and we find it to be too small to 
produce any observable effect on the Sun properties. 

However, for scenarios with null or very small anni- 
hilation cross sections, (crv) ;$ 10 -33 cm 3 s _1 , the num- 
ber of trapped WIMPs is significantly increased. In this 
case, the transport of energy by the WIMPs from the 
interior of the Sun to the outer shells can dramatically 
reduce the core temperature, the most important conse- 
quence being a reduction of the S B flux. Considering the 
present theoretical and experimental uncertainties on §" B , 
we have studied the combination of DM masses and SD 
scattering cross sections which can be ruled our with this 
argument, finding only a small region of the parameter 
space which is not already excluded by direct detection 
bounds. Even with a significant decrease of the uncer- 
tainties on <j) v B , the region of the parameter space which 
can be probed remains approximately the same, there- 
fore future experimental advances will not significantly 
change the situation. 

We have also considered SIDM candidates which have 
been recently invoked as a solution of the solar compo- 
sition problem. Correctly implementing the geometrical 
limit, we find that the number of WIMPs captured on the 
star is not sufficient to produce significant deviations on 
the Sun properties, unless one considers very high scat- 
tering cross sections, where however the models become 
in tension with direct detection bounds. Even in this 
case, we encounter significant modifications of the den- 
sity and temperature profiles only in the inner regions 
of the Sun while the outer shells arc not significantly af- 
fected. We can therefore exclude that the transport of 



energy produced by a population of WIMPs can solve or 
even alleviate the solar composition problem. 

Finally, we notice that the WIMPs energy transport 
may produce dramatic effects on the structure of stars 
settled in high DM density environments, like the galactic 
center. We leave the investigation of this scenario for 
future work. 

When this manuscript was ready for submission, a sim- 
ilar work appeared on the arXiv, focusing on the effects 
of WIMPs in the Sun [63. While our work mainly fo- 
cus on the modifications of the solar neutrino fluxes, Ref. 
J69t | considers future helioseismology data as a diagnostic 
tool to constrain the WIMPs parameter space. 
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Appendix A: Formalism and code 

1. Capture 

The formalism to compute the rate at which WIMP 
particles are captured by a star have been extensively 
studied in the eighties. Here we adopt the results from 
27 1, who gives for the capture rate C: 



C = V] 47r / drr 



.dCj{r) 
dV 



(Al) 



with 



dC t {r) /6 



dV 



1/2 



Pi(r) Px v2 ( r ) 



a *' Ni ' Mi m x v 2 2r)A 2 



(A2) 



\[A + A---)[ X (-V,V)-X(A-,A + )] 



1 a -A 2 1 , -A 2 1 -n 2 

-A+e - - -A_e A + - -Tie n 

2 + 2 2 ' 



2 _ 3^ 2 (r) M 



A z = 



2v 2 p 2 _ ' 



A± = A ± 7), rf 



3^ 
2v 2 



11 



X(a, b) = ^[Erf(fe) - Erf(o)] = / dye 



\x- = {pn - l)/2, p t = m x /Mi 

where pi (r) is the density profile of a given chemical ele- 
ment in the interior of the star and Mi refers to its atomic 
mass, while m x and p x are respectively the WIMP mass 
and the WIMP density at the star position. The analytic 
expression for the capture rate per shell volume reported 
above, dC/dV, is obtained for the case of a Maxwell- 
Boltzmann velocity distribution with speed dispersion v; 
throughout this paper, we adopt v = 270 km s~\ and 
we take the velocity of the star moving through the DM 
halo, labeled as v*, to be v* = 220km s^ 1 , as appro- 
priate for the Sun. The radial escape velocity profile 
depends on M(r), i.e. the mass enclosed within a ra- 
dius r, v 2 (r) = 2j™GM(r')/r' 2 dr'. The DM elastic 
scattering cross section off nuclei cr X Ni is the sum of the 
Spin-Dependent (SD) and Spin-Independent (SI) contri- 
butions: 



T xNi 



f3 z Af<j SI + P'osd 



4(J< + 1) 
3Ji 



\(Si, P ) + (S hn )f 



with osi and osd the hydrogen-normalized SI and SD 
nuclear-scattering cross sections. The factor Aj is the 
atomic number while J, is its spin and j3 the ratio of the 
reduced mass of the WIMP-nucleus and WIMP-proton 
systems. Finally, (5j iP ) and (Si >n ) are respectively the 
expectation values of the spin content of the proton and 
neutron group in the nucleus. The total capture rate is 
then obtained summing up all the elements present in 
the star. 

The geometrical size of the star itself sets an up- 
per bound to the effective WIMP scattering cross- 
section, <J & g, which we model imposing a maximum to 
the capture rate once the condition below is reached: 



T eff 



E° 



xNi 



Ni = nRi 



with Ni the number of nuclei of the element i; this is 
sometime referred to as the optical limit. 

2. Annihilation 

Throughout this paper we work under the assump- 
tion that WIMPs thcrmalizc inside the Sun within a 
negligigible timespan with respect to stellar evolution- 
ary timcscalc. This is a good approximation for most 
WIMP models, and especially for the high scattering 
cross-sections we will examine, as an upper limit of the 
thermalization time reads 11311: 



4-7T 



Tth 



R 7/2 
-ft* 



>2G v xNi M, 



3/2 



(A3) 



that for the sun means 



Tth 



10 5 s 



10" 36 cm 2 



lOGcV 



a xNi 



In the case of a non-local WIMP energy transport, i.e. 
for large K (see next section), the WIMPs clould shares 
a global temperature which can be set, in good approx- 
imation, to the star core temperature T c . The WIMP 
radial distribution in the star gravitational well is thus 
given by [79| : 



n x ,nl(r) = n xfi e~^~ ,r x 



3kT r 



2TrGp c m x 



(A4) 



with p c referring to the core density. The distribution 
results quite concentrated toward the center of the star: 
typical values for ZAMS stars are r x ~0.1 r c ~ 10 9 cm, 
r c being the stellar core radius. Such zero order distribu- 
tion is slightly modified by transport effects, see Section 
I A 41 and we adopt the modified distribution to compute 
the actual rate of annihilations; however, the effects of 
such modification are small and their effects on annihi- 
lations negligible in practice. Once the normalization 
n Xt o is obtained solving the Eqjl] the distribution n x (r) 
is completely specified and the annihilation term can be 
easily computed: 



/•Re 
Jo 



A = / e ann r <±np(r)dr, 



with: 



e<mn = 7;\< JV ) rn x c n x 



(r)p(r 



(A5) 



(A6) 



p(r) is the baryon density at a given position and e ann 
is the luminosity produced by WIMPs annihilations per 
unit of baryonic mass. The factor 1/2 (1/4) in the equa- 
tion above is appropriate for self (not self) conjugate par- 
ticles and (<jv) the velocity-averaged annihilation cross- 
section. If an equilibrium between capture and annihila- 
tion is reached the annihilation rates reduces to A=C/2 
and it is independent on the annihilation cross-section. 
The equilibrium timescale for such a processes, neglect- 
ing evaporation, is given by 



(J- 
\CA 



1/2 



V^r 3 / 
C(av) 



1/2 



(A7) 



and for annihilations cross-sections of the order (av) ~ 
10 -26 cm 3 s _1 , a typical value suggested by the require- 
ment of the correct WIMP relic density, it is much shorter 
than the age of the Sun. In our calculations we compute 
the capture and annihilation terms self-consistently, and 
find good agreement with the equilibrium approximation 
for this value of the annihilations cross-section. However, 
we also explore scenarios with negligible or null annihila- 
tion cross-sections, as for the case of asymmetric- WIMPs, 
for which the equilibrium between annihilations and cap- 
ture is not reached. 
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3. Evaporation 

WIMPs trapped inside the Sun may scatter off nuclei 
to velocities high enough to escape the gravitational field 
of the star and thus leave the system. Detailed studies on 
this phenomena [80, HlJ , called evaporation, have shown 
that the evaporation timescale, defined as the inverse of 
the evaporation rate, E in Eq[T] depends exponentially 
on the mass of the WIMPs. Therefore, below a certain 
mass threshold, m ev , basically all the WIMPs evaporate 
during the Sun evolution and conversely, slightly above 
m ev the evaporation is negligible. At high WIMP scatter- 
ing cross sections, corresponding to the LTE regime, the 
evaporation starts to be inefficient since the WIMPs leav- 
ing the center of the star are rapidly re-scattered to low 
velocity orbits. The evaporation mass, depends therefore 
of the size and nature (SI or SD) of the scattering cross 
section. Considering the Sun, for agD ~ 10~ 36 cm 2 it 
is typically of the order m ev ~ 4 GeV and it decreases 
to m ev ~ 2 GeV at <tsd ~ 10~ 34 cm 2 . In the following 
we focus to WIMP masses m x > 5 GeV so the effect 
of the evaporation can be completely neglected for our 
purposes. 



x exp 



ka(r 



dr' 



d<p(r') 



"X dr 



kT*(r') 



— dr' 



with the normalization n x ,LTE(0,i) obtained imposing 

J * 47rr 2 n XiLTE (r, t) dr = N(t). T r (r) is the tempera- 
ture profile of the star and </>(r) is the local gravitational 
potential. The thermal diffusivity and conductivity coef- 
ficients, respectively n(r) and a(r), are defined as: 



i ^3 a 3 



(JiTiiir) 



nj{r,t) 



ai{iii) 



and 



<r) = {l(r)Y,H»iMr)} '} \ 



(A12) 



(A13) 



with OjXni) and Kj(jUj) obtained from tabulated values in 
Ref. |82|. The sum runs over all the nuclear species, Hi is 
the WIMP-to- nucleus mass ratio and n., (r) is the number 
density of the species i at any given radius. 

The two correction factors which appear in Eq. IA8I 
read: 



4. Transport 

The energy (per unit mass, per unit time) transported 
by WIMPs can be cast as: 



1 



etr 



with 



47rr 2 p(r) dr 



[WMr)L 



trans, LTE 



(r)]; (A8) 



Lti £ 



,LTE 



(r) =47rr 2 K(r)n x x TE (r)Z(r) 
kT*(r)-|V2 k _di;(r) 



(A9) 



dr 



being the "transport luminosity" to be imputed to the 
WIMPs inside the Sun. The Knudsen number, K, indi- 
cates the "localization" of the WIMPs transport: 



K = 



1(0) 



(A10) 



The WIMP distribution ?i X: lte i s obtained in the ap- 
proximation that the Knudsen number is very small and 
WIMPs are in local equilibrium with the baryons: 



™ x ,LTE(r) = 7-i XiLT e(0) 



r*(r) 



T, 



3/2 



(All) 



H^c-^y + 



(A14) 



and 



KK 



l- 



i 



-{\nK-\nK )/i 



= 1- 



1 



. K 



,iA 



1-1- e v""- '"- "i> ■ 

' (A15) 
with t = 0.5 and K o =0A. 

These two quantities, introduced by Gould & Raffelt 
[82| , extend the treatment of the WIMPs energy trans- 
port obtained under the assumption of local thermal 
equilibrium (i.e. small K), to the case of large Knud- 
sen numbers. Therefore, the equations here described 
are virtually correct for any value of WIMP scattering 
cross-sections, applying both to local and non-local en- 
ergy transport regimes. The WIMPs energy distribution 
in Eq. IA11I differs from the analogue expression Eq. IA4I 
valid for large K. In the following, for the computation of 
the annihilation rate, we adopt a distribution interpolat- 
ing between the two regimes that we defined, following 
0, as: 

n x( r ) = /(-K>x,LTK(r) + (1 - f(K))n x , n i. 



As commented in Sec lA2l this modification is however 
negligible in practice. 
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